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CDF Collaboration 
Abstract 

We present a study of the production of K° s and A in inelastic pp collisions at \fs= 1800 and 630 GeV 
using data collected by the CDF experiment at the Fermilab Tevatron. Analyses of K® and A multiplicity 
and transverse momentum distributions, as well as of the dependencies of the average number and (pt) of 
Kg and A on charged particle multiplicity are reported. Systematic comparisons are performed for the 
full sample of inelastic collisions, and for the low and high momentum transfer subsamples, at the two 
energies. The pr distributions extend above 8 GeV/ c, showing a {px) higher than previous measurements. 
The dependence of the mean Kg (A ) px on the charged particle multiplicity for the three samples shows 
a behavior analogous to that of charged primary tracks. 



PACS numbers: 13.85.Hd, 13.87.Fh 



I. INTRODUCTION 

Hadron interactions are often classified as either "hard" or "soft" . Although there is no 
formal definition for either, the term "hard interactions" denotes high momentum transfer parton- 
parton interactions typically associated with such phenomena as jets of high energy transverse to 
the incoming hadron momenta (Ex). The "soft" interaction component encompasses everything 
else and dominates the inelastic cross-section. From a theoretical point of view, perturbative QCD 
provides a reasonable description of high-.&r jet production. However, non-perturbative QCD, 
relevant to Iow-Et hadronic production, is not well understood. Some QCD-inspired models ^ 
attempt to describe these processes by the superposition of many parton interactions extrapolated 
to very low momentum transfers. It is not known, however, if these or other collective multi-parton 
processes are at work. The experimental studies of low-.&r interactions are usually performed on 
data collected using minimum bias (MB) triggers, which, ideally, sample events in fixed proportion 
to the production rate — in other words, in their "natural" distribution. Lacking a comprehensive 
description of the microscopic processes [3( involved in low-.&r interactions, our knowledge of the de- 
tails of low transverse momentum (p?) particle production rests largely upon empirical connections 
between phenomenological models and data collected with MB triggers at many center-of-mass en- 
ergies (E cms ). Such comparisons necessarily face the difficulty of isolating events of a purely "soft" 
or "hard" nature. 

Comparative studies of the event structure through collective variables such as the charged 
particle multiplicity and the transverse energy of the event are important to our understanding 
of the soft production mechanism. In a previous paper 4], a novel approach in addressing this 
issue using samples of pp collisions at = 1800 and 630 GeV collected with an MB trigger was 
described. The analysis divided the full MB samples into two subsamples, one highly enriched in 
soft interactions, the other in hard interactions. Comparisons between the subsamples and as a 
function of E cms were performed. The same approach has been applied here to the production of 
strange particles. 

Beside gluons and the lighter quarks u and d, strange quark production is the only component 
of \ow~pt multiparticle interactions which is statistically significant and experimentally accessible 
with an MB trigger. It is also a probe for investigating the transition of soft hadron interactions to 
the QCD high-p-r perturbative region. 

This paper describes a study of and A production in pp interactions at different E cms . 



Inclusive distributions of the multiplicity and transverse momentum of K® and A are presented 
first. The high statistics of the data sample collected at \^s= 1800 and 630 GeV allow an extension 
of the range and precision of these measurements with respect to previous ones. Studies of the 
dependence of the average pt of iff (A ) and of their mean number on the event charged multiplicity 
are also presented. Different behavior of the hard and soft subsamples is observed, consistent with 
prior reports on charged particles Q. 



II. DATA COLLECTION 
A. The CDF Detector 



Data samples have been collected with the CDF detector at the Fermilab Tevatron Collider. 

n 

The CDF apparatus has been described elsewhere [5J; here only the parts of the detector utilized 
for the present analysis are discussed. The coordinate system is defined with respect to the proton 
beam direction, which defines the positive z direction, while the azimuthal angle <fi is measured 
around the beam axis. The polar angle 9 is measured with respect to the positive z direction. The 
pseudorapidity, rj, is often used and is defined as rj = — ln(tan[#/2]). Transverse components of 
particle energy and momentum are conventionally defined as projections onto the plane transverse 
to the beam line, = E sin 6 and pt = \p\ sin 9. 

Data were collected with an MB trigger at 1800 GeV during Runs 1A (1992-93) and IB (1994- 
95), and at 1800 and 630 GeV during Run 1C (1995-96). This trigger requires coincident hits in 
scintillator counters, located at 5.8 m from either side down stream of the nominal interaction point 
and covering the pseudorapidity interval 3.2 < \q\ < 5.9, in coincidence with a beam-crossing. 

The analysis uses charged tracks reconstructed within the Central Tracking Chamber (CTC). The 
CTC is a cylindrical drift chamber covering an 77 interval of about three units with high efficiency for 
\q\ < 1 and px > 0.4 GeV/c. The inner radius of the CTC is 31.0 cm and the outer radius is 132.5 
cm. The full CTC volume is contained in a superconducting solenoidal magnet which operates at 
1.4 T The CTC has 84 sampling wire layers, organized in 5 axial and 4 stereo "superlayers" Q|. 
Axial superlayers have 12 radially separated layers of sense wires, parallel to the z-axis (the beam 
axis), that measure the r-0 position of a track. Stereo superlayers have 6 sense wire layers, with a 
~ 3° stereo angle, which measure a combination of r-<p and z positions. The stereo angle direction 
alternates with each neighboring stereo superlayer. Measurements from axial and stereo superlayers 



are combined to form a three-dimensional track. The spatial resolution of each point measurement 
in the CTC is less than 200 /jm; the transverse momentum resolution, including multiple-scattering 
effects, is a PT /p 2 T < 0.003 (GeV/cT 1 . 

Inside the CTC inner radius, a set of time-projection chambers (VTX) 8] provides r-z tracking 
information out to a radius of 22 cm for \rj\ < 3.25. The VTX is used in this analysis to find the z 
positions of event vertices, defined as sets of tracks converging to the same point along the z-axis. 
The closest detector to the beam-pipe is the Silicon Vertex Finder (SVX), used to reconstruct ver- 
tex positions in the transverse view. Reconstructed vertices are classified as either "primary" or 
"secondary" based upon several parameters: a minimum of 4 converging track segments in \rj\ < 3 
(a track segment is a sequence of 4 aligned hits), the total number of hits used to form a seg- 
ment, forward-backward symmetry and vertex isolation. Isolated, higher multiplicity vertices with 
highly symmetric topologies are typically classified as primary; lower multiplicity, highly asymmet- 
ric vertices, or those with few hits in the reconstructed tracks, are typically classified as secondary. 
Systematic uncertainties introduced by the vertex classification scheme are discussed in Section IVT1 

The transverse energy flux is measured by a calorimeter system [9] covering \rj\ < 4.2. The 
calorimeter consists of three sub-systems, each with separate electromagnetic and hadronic compo- 
nents: the central calorimeter, covering the range \rj\ <1.1; the end-plug, covering 1.1< \r)\ <2.4; 
and the forward calorimeter, covering 2.2< \r)\ <4.2. Energy measurements are made within pro- 
jective "towers" that span 0.1 units of rj and 15°(5°) in within the central (end-plug and forward) 
calorimeter. 

B. The Data Set 

The 1800 GeV MB data sample consists of subsamples collected during three different time 
periods. Approximately 1.7xl0 6 events were collected in Run 1A at an average luminosity of 
3.3xl0 30 s _1 cm~ 2 , 1.5xl0 6 in Run IB at an average luminosity of 9.1xl0 30 s _1 cm~ 2 and 1.06xl0 5 
in Run 1C at an average luminosity of 9.0xl0 30 s _1 cm~ 2 . The 630 GeV data set consists of about 
2.6xl0 6 events recorded during Run 1C at an average luminosity of 1.3xl0 30 s^'cm" 2 . 

Additional event selection conducted offline removed the following events: (i) events identified 
as containing cosmic ray particles as determined by time-of-flight measurements using scintillator 
counters in the central calorimeter; (ii) events with no reconstructed tracks; (iii) events exhibiting 
symptoms of known calorimeter problems; (iv) events with at least one charged particle recon- 



structed in the CTC to have pt > 400 MeV/c, but no central calorimeter tower with energy 
deposition above 100 MeV; (v) events with more than one primary vertex; (vi) events with a pri- 
mary vertex more than 60 cm away from the center of the detector (in order to ensure uniform 
acceptance in the assumed fiducial region and good track and calorimeter energy reconstruction); 
and (vii) events with no primary vertices. 

After all event selection requirements, 2,079,558 events remain in the full MB sample at 
y/s = 1800 GeV and 1,963,157 at yfs = 630 GeV. The vast majority of rejected events failed 
the vertex selection. About 0.01% of selected events contain background tracks from cosmic rays 
that are coincident in time with the beam crossing and pass near the event vertex. The residual 
contamination due to the interactions of the beam particles with the gas in the beam pipe is about 
0.02%. A more detailed discussion of the systematic uncertainties arising from the event selection 
criteria and other sources is presented in Ref. Q|. 



III. CHARGED TRACKS AND i^(A°) SELECTION 

We require all reconstructed tracks to pass through a minimum number of layers in the CTC and 
have a minimum number of hits in each superlayer in order to reduce the number of misreconstructed 
tracks and those with large reconstruction uncertainties. The remaining track set, which includes 
primary and secondary tracks, is used as a starting point for both the selection of primary charged 
tracks and for the A, (A ) candidate identification procedure. 

Charged track multiplicity definition. Tracks are required to pass within 0.5 cm of the 
beam axis, and within 5 cm along the z-axis from the primary event vertex. In order to ensure high 
efficiency and acceptance, tracks are accepted only if they satisfy the conditions pr> 0.4 GeV/c 
and \tj\ < 1.0 . This selection defines the charged track multiplicity in an event, N* h . 

A s ° and A selection. A s ° and A jiol] ffrom now on collectively referred to as V^ ) are selected 
looking for opposite-charge pairs of tracks converging to a common vertex displaced from the beam 
line in the transverse direction. A vertex fit is performed to ensure that the two tracks originate 
from the same vertex. A candidate is required to have a fit probability greater than 5%. In a further 
step a fit is performed constraining the V° momentum vector (within the track uncertainties) to 
point in the direction of the primary vertex (pointing constraint fit). The candidates are kept if 
the fit probability is greater than 5% and the recomputed invariant mass is within three standard 



deviations of the word average K® or A mass 
The analysis selection also requires: 

• L xy (V°)>l cm, where L xy is the distance from primary vertex to the decay vertex of the V° 
in the r-0 plane; 

• both decay tracks have \r)\ < 1.5 and pt > 0.3 GeV/ c; 

• the V° line-of- flight is close to the event vertex along the z axis: \zq° — z^ ertex \ <6 cm; 

• impact parameter do(V°) < 0.7 cm; 

• Pt(V°) >0.4 GeV/cand \r](V°)\ <1.0 

For events with more than one V° candidate sharing the same track, only the candidate with 
the lower vertex fit x 2 /Nd.oj. is retained. 

After all selection requirements, we find 36,642 and 7,518 A in the 1800 GeV MB sample 
and 32,222 and 5,883 A in the 630 GeV MB sample (see Table [Q). 

The invariant mass distributions of the and A surviving the selection requirements, but with 
the mass window extended to ten standard deviations from the world average, are shown in Fig. ^ 
in both cases the peaks are narrow but, because of the fit procedure, the background is not flat and 
may not be accounted for by the level of the sidebands. We also note that this background includes 
the contamination of K® in the A sample and vice versa. A detailed background evaluation is 
discussed in Section IVl 



IV. SELECTION OF SOFT AND HARD INTERACTIONS 



12 1 since 



The identification of "soft" and "hard" interactions is largely a matter of definition 
it is unknown how to distinguish soft and hard parton interactions. This is true from both the 
theoretical and experimental points of view. In this analysis, we use a jet reconstruction algorithm 
to define the two cases. The algorithm employs a cone with radius R = (At/ 2 + A0 2 ) 1/2 = 0.7 to 
define "clusters" of calorimeter towers belonging to a jet. To be considered, a cluster must have a 
transverse energy Et, defined as the scalar sum of the transverse energy of all the towers included 
in the cone, of at least 1 GeV in a seed tower, plus at least 0.1 GeV in an adjacent tower. 

In the regions |r/| < 0.02 and 1.1 < |?7| < 1.2, a track-clustering algorithm is used instead of the 
calorimeter algorithm to compensate for energy lost in calorimeter cracks. A track cluster is defined 



as one track with > 0.7 GeV/ c and at least one other track with pt >0.4 GeV/ c in a cone of 
radius R = 0.7. 

We define a soft event as one that contains no cluster with Et > 1.1 GeV. All other events are 
classified as hard. 

V. EFFICIENCY AND CORRECTIONS 

The probability of observing a real V° in the apparatus is influenced by several effects. In this 
section we discuss the efficiency of track reconstruction, the correction for limited acceptance, and 
evaluation of the background. At the end, some cross-checks of the correction procedures are also 
briefly described. 

1. The efficiency for finding Kg(A°) has been investigated in two different ways. In the first 
method, simulated hits from singly-generated V° are embedded among the set of hits of MB events 
from the data. The events are then reconstructed with standard V° search and selection. In the 
second method, entire MB events with V° production and decay are generated with pythia/jetset 
Monte Carlo (MC) j^J, Q|. Full CDF detector simulation and reconstruction are then applied to 
the events and the resulting reconstructed kinematic distributions are similar to those observed in 
the data. The results from the two methods are compatible within the statistical uncertainties. 

The efficiency is defined as the ratio of the reconstructed to the generated number of V° in the 
fiducial region. It is examined as a function of single kinematical variables of the V°, integrating 
over all the remaining variables. The embedding method, given its almost flat V° distribution 
in all variables, gives smooth and statistically better determined efficiency dependences from all 
observables over all the acceptance limits. The results of this method are used to determine the 
shape of the efficiency as a function of any chosen variable. Each efficiency distribution from the 
embedding method is then scaled by an overall normalization factor so that the integrated efficiency 
obtained from the embedding method matches the integrated efficiency from the full MC method. 

The efficiency for finding a single K®(A°) is approximately constant (around 40% (32%)) as a 
function of r](V°) in the region of |?y(V°)| <1 and pt{V°)>0A GeV/c. As a function of pr{V°) (in 
the same r\ region), the efficiency rises rapidly from 25% (15%) at 0.4 GeV/c to about 50% (40%) 
for p T ~l GeV/ c, and then slowly decreases to ~ 20% (~ 15%) for p T > 8 GeV/ c. This behavior is 
due to the difficulty in reconstructing low-pr secondary tracks and in identifying secondary vertices 
far from the primary vertex. The efficiency also diminishes for L xy < 3 cm, while it is roughly 



constant as a function of the charged multiplicity of the event. The overall efficiency is about 39% 
for and 31% for A . 

2. A correction for the fiducial acceptance requirement in L xy and in the pr of the K®(A°) decay 
products is estimated using MC and found to range from about 15(20) at pr= 0.5(1) GeV/c to 
about 1 for p T > 5 GeV/ c. 

3. The contamination by A in the sample is estimated to be ~ 3% as found in the pythia 
MC simulation; the contamination by K® in the A sample is about 7% on average while it is almost 
50% for p T (A°) < 1.5 GeV/ c. The same MC sample is used to compute the probability of selecting 
fake secondary vertices (not due to or A decays). Such probability is found to account for 
roughly 25% of the and 40% of the A . 

4. The overall correction factor for a generic inclusive variable X (e.g. the Pt{V )) is given by 
the expression: 

C(X) 



1 — Rfake(X) 



e(X) x A(X) 

where Rf a ke is the probability of a fake V°, e is the global efficiency and A is the acceptance. 
The overall correction factors, as a function of Pt{V°), are shown in Fig. El The integrated MC 
correction factors are estimated to be 4.5±0.1 and 10.1±0.2 for K® and A respectively. 

5. Because of the small differences that exist between some pythia distributions and the data, 
we expect that the MC correction will not be fully reliable in the regions where it changes very 
rapidly. Evidence of this is given by the reconstructed V° px versus the proper time which shows 
a depletion in the low-pr and low-lifetime region, even after applying the MC correction. 

We use the following method to correct the counted number of i^°(A ) in this region. The K^{hP) 
invariant px distributions for the full MB sample are fitted with a functional power-law form: 

dp 6 \po + p T J 

where E is the particle energy and po, A and n are free parameters, in the region above 0.8 GeV/ c 
(1.1 GeV/ c for A ). This equation has been widely used to fit the pt distributions of charged tracks 
down to the lower measured pt [l^. The fitted function is extrapolated down to pt(V°)=0A GeV/ c 
and the corrected number of K®(A°) is extracted from the integral of the curve. 
In the full MB sample, the number of undetected V° is estimated to be approximately 18 x 10 3 K® 
and 14 x 10 3 A at 1800 GeV, and 24 x 10 3 and 12 x 10 3 , respectively, at 630 GeV. 

6. The above correction affects the measurements of the mean number of V° per event and 
of the mean p T when computed at fixed N* h . The latter is calculated as the sum of the pr's of 



Kg (A ), above 0.4 GeV/ c, observed in events of a given charged multiplicity, divided by the number 
of ir s °(A°): 

I N vo 

(Pt) = -Tj— Pr t (3) 

An estimate of the number of undetected V° and the resulting effect on the (pr(^ )) and on the 
V° multiplicity are obtained, for each pt{V°) distribution, with the procedure used in the inclusive 
case. The constraint that the sum of undetected V° for each multiplicity should give the number 
of undetected V° computed from the inclusive pt distribution is imposed. The corrected (pt{V )) 
is computed by extrapolating the fitted Pt{V°) distribution down to pr(^ )=0.4 GeV/ c. 

7. The consistency of the correction procedures described above has been verified through the 
following cross-checks. 

In order to check the selection requirements and the quality of the efficiency correction, the 
raw and corrected proper lifetime distributions at 1800 GeV are shown in Fig. El Fitting to an 
exponential form gives a mean proper lifetime of (0.89 ±0.01) x 10~ 10 s (x 2 /^d.o./.=49.7/59) for 
1800 GeV and (0.90 ±0.01) x 10" 10 s (x 2 /N d . oJ =60.6/56) for 630 GeV. Both values are consistent 
with the world average values ^l|. The same fit to the A proper lifetime distributions gives a 
mean of (2.61 ±0.07) x 10~ 10 s {x 2 /N dM . f =44.2/49) for 1800 GeV and (2.61 ±0.07) x 10~ 10 s 
(x 2 /A fio ./.=57.4/50) for 630 GeV. The proper lifetime regions used for the fit are r > 0.7 x 10~ 10 s 
(K°) and t > 10- 10 s (A ). 

The number of undetected V° extracted from the fitted pt curve is also checked. The proper 
lifetime distributions of Fig. El are fitted to an exponential form with fixed slope (the K®/A° mean 
lifetimes Q, t k o = 0.8935 x 10~ 10 s; r A o = 2.632 x 10~ 10 s) in the region r > 0.7 x 10~ 10 s 
(r > 10 -10 s for A ) and the fitted curves are integrated down to r = 0. The number of undetected 
K®(A°) obtained matches to 15% (30%) with the number from the pt distribution. Furthermore 
the pt distributions of K®(A°) with proper lifetimes greater than 0.8 x 10~ 10 s (1.0 x 10~ 10 s for A ) 
are compared with the corresponding distributions for all lifetimes; the comparison gives the same 
values of average pt- When normalized to one another, the curves give a comparable number of 

(A ) in the extrapolated region. 

An additional cross-check for correcting the average pt of the V° observed in events of fixed mul- 
tiplicity consists of plotting the proper lifetime distribution in slices of pr so that each distribution 
corresponds to one bin in p?. This is done for each bin in multiplicity. After fitting the distribution 
in the long lifetime region in each pt bin, the correct number of K® in the short lifetime region 



can be extrapolated from the fits. The (pt{V )) can then be recomputed from the modified pr 
distribution. The px values obtained using the two different correction methods are consistent. In 
the case of A , no events are found with px below 1 GeV/ c due to the tight fiducial requirements 
imposed in the analysis. Therefore, in the A case, the correction method based on extrapolating 
the proper lifetime distribution at each p? bin cannot be used. Because of this, the cross-checks are 
limited to comparing the number of extrapolated A in the px and proper lifetime distributions of 
the full data sample. 

n 

Finally, we refer to [4| for a detailed discussion of the charged track selection and reconstruction 
efficiencies. 

VI. SYSTEMATIC UNCERTAINTIES 

The two dominant systematic uncertainties come from the acceptance and efficiency correction 
procedures. As described in Section El acceptance and efficiency corrections have been computed 
using MC simulation, with an additional correction applied to compensate for MC deficiencies in 
the low pt region between 0.4-0.8 GeV/ c. The two correction procedures are largely independent, 
which allows us to evaluate the systematic uncertainties from these two sources separately. The 
details are described below. 

1. We study the sensitivity of this measurement to the differences between the MC predictions 
and the shapes of the observed V° kinematical distributions. We use the following two sets of MC 
events. The first is created using the default pythia MC. The second is the one used for efficiency 
studies using the embedding procedure: the V°'s in this set have non-physical distributions roughly 
uniform in p? but not in 7]. The different correction factors evaluated from the above data sets are 
applied to the measured distributions. Half the difference between the corrected distributions is 
taken as the systematic uncertainty on the distributions themselves, which amounts to about 10% 
for the and A px distributions, roughly constant over the whole spectrum. 

The effect on the mean pt value is 3% for and 4% for A . For the and A multiplicity 
distributions, the systematic variation ranges from 10% to ~25%. As a function of N* h , the sys- 
tematic uncertainty on the number of ranges from a few percent to roughly 20% at the highest 
charged multiplicities. 

2. The systematic uncertainty due to the correction for the undetected V° in the pt region 



between 0.4 and 0.8 GeV/ c has been evaluated in the following way. The procedure defined in 
Section |V[ point 7, is repeated using pt and proper lifetime distributions both corrected with 
pythia MC and with the embedding-based correction. The total number of V° is computed by 
integrating the corrected pr and lifetime spectra for each of the two cases. We end up with four 
different evaluations of the number of undetected V°. By comparing the numbers obtained from 
all combinations, we observe that the largest difference amounts to about 50% of the correction 
value. This number is taken as the systematic uncertainty on this correction and is counted as a 
contribution to the systematics on the the total number of measured if s °(A°). 

The mean p? values at fixed multiplicities are also affected by the correction for the undetected 
if,? (A ) in the low pt region. The systematic uncertainty on the correction is estimated as follows. 
First it has been verified that the mean p? after correction is independent of the K®(A°) proper 
lifetime in the region used in this analysis (see Fig.0J). Then, starting with p? distributions at fixed 
charged multiplicity for the subset of events with K®(A°) proper lifetime greater than 0.8 x 10~ 10 s 
(1.0 x 10~ 10 s), the mean p? is computed the same way as described in Section IS and the difference 
between the mean p T values for the full dataset and the high r subset is assigned as a systematic 
uncertainty for this correction. Since the correction is applied only to calculations of the mean 
Pt at fixed multiplicity and of the number of iff (A ), the systematic uncertainty associated with 
it affects only these measurements. It amounts to about 6% (10%) for the total number of A s ° 
(A ) and affects the average number of K®(A°) as a function of the charged multiplicity by the 
same amount. These systematic uncertainties combined in quadrature with the other systematic 
uncertainties discussed in this section are included in Figures [TT] to UHl 

3. To investigate the systematic effect of the track reconstruction procedure on the efficiency 
correction, we compare our result with a set of MC events where the tracks are reconstructed using 
the CTC information alone, as opposed to the default SVX+CTC track reconstruction. We find 
that the variation on the final corrected pt distribution is negligible. 

4. Other sources of systematic uncertainties include the dependence of the results on the in- 
stantaneous luminosity and the uncertainty associated with the identification and selection of good 
isolated pp interactions from secondary or closely spaced event vertices (see j^j]). The first may 
affect the results because higher luminosity gives higher detector occupancy which in turn can alter 
the V° identification. This has been investigated by analyzing data samples recorded at different 
instantaneous luminosities. The results show no observable effect. The second source can lead to 



incorrect event selection and produce associations of tracks that fake a V . This source has been 
investigated by comparing data samples with different requirements for a good pp vertex 0] . The 
results give systematic variations smaller than 9% on the overall number of K®(A°). 

5. The uncertainty on the correction is distributed in different ways for different observables. As 
a consequence, the integral of the corrected distribution of each variable is different. For example, 
the total number of iff (A ) extracted from the integral of the corrected pt distribution may be very 
different from that extracted from the multiplicity distribution. In particular, as discussed in the 
previous section, the px correction has been observed to be unreliable for p? < 0.7 GeV/c where 
a large part of the V° cross-section lies, so that the area under the distribution may be subject 
to large uncertainties. Given this, we use the global (integrated) correction from the pythia MC 
as a correction factor for the total number of K®(A°). We renormalize each distribution to this 
number to which we attribute a 30% systematic uncertainty. This value is determined as the 
maximum difference that was found between the global corrected number of V° and the integral of 
any corrected distribution. Such uncertainty reflects on the K^/ir ratios and on the absolute scale 
of the ratios of the mean number of A^°(A°) to the charged multiplicity plotted in Figures IT7I to 1201 
Table ITTT1 reports a summary of all the systematic uncertainties discussed. 

VII. ANALYSIS RESULTS 
A. Results 

All data presented are subject to > 0.4 GeV/c and \r]\ < 1 requirements, as specified in 
Section UTT] and are corrected for acceptance and vertex-finding efficiency. Systematic uncertainties 
are not included except where explicitly stated. Table |l] shows the raw and corrected numbers of 
A^°(A°) selected in our fiducial region for the full MB sample as well as for the soft and the hard 
samples. The corrected mean number of K®(A°) per event in each sample is also shown; systematic 
uncertainties are included. 

In Fig. for the K® and in Fig. [B] for the A , the normalized multiplicity of K®(A°) for the 
MB, soft and hard events is shown separately for the ^=1800 GeV (solid symbols) and 630 GeV 
(open symbols) data. The probability of producing one or more A is lower than the equivalent 
probability, and the difference increases with V° multiplicity. This behavior is more pronounced in 
the soft subsample. The results shown in Figures 03 and El with their statistical errors, are reported 



in Table 

The invariant pt inclusive distributions of K® are shown in Figures [7|and|Hlat the two energies 
for the full MB, soft and hard samples. Data are normalized to the number of events in each 
sample. Figures M and ITU1 show the same px distributions for the A . 

The dependence of the K® and A average px, calculated as described in Eq. (jSJ), on the event 
charged multiplicity is shown in Figures fTTl- T 1800 GeV) and El-El (630 GeV). The mean p? of 
primary charged tracks measured in the same phase space region, as published in is also shown 
for comparison. For the dataset, in the region ranging from 0.4 to 0.8 GeV/ c, the corrected data 
points are assumed to lay on a curve of form (J2J) extrapolated from the fit to the measured data 
points in the region px >0.8 GeV/ c (details of the correction procedure are described in Section IV)). 
Note, that with the kinematical selection used in this analysis, no events with p^(A°)<l GeV/c 
were observed. For the measurement of (pr), we fit the spectrum in the region of pt>1.1 GeV/c 
using Eq.© and extrapolate down to pr=0.4 GeV/ c. We define the (pt) as the mean value of the 
fitted function. This definition is adopted in order to compare the (pt) with that of and of 
charged tracks. 

Figures El to [201 show the ratio of the mean number of K®(A ) per event to the multiplicity as a 
function of the multiplicity itself. The charged particle multiplicity N* h was chosen as the reference 
variable to analyze V° production. The reason for this choice is based on the observation that 
the event charged multiplicity is a global event variable characterizing the whole multiparticle final 
state and is related to the hardness of the interaction (see [3). As in the case of charged 

particles, possible new structures in the V° final state correlations would be exhibited as a function 
of N* h . The dependence of the average px on multiplicity, for example, remains unexplained in any 
of the current models. 



B. Dependence on E? threshold 

It has been remarked in the previous sections that the identification of soft and hard events is 
essentially a matter of definition. In order to investigate the sensitivity of the above results to the 
the cluster energy threshold used to separate soft and hard events, the analysis has been repeated 
changing the threshold from 1.1 to 3.0 GeV. Although, as expected, the higher threshold value 
influences the global statistics of the soft and hard components, it preserves the shapes of the 
inclusive pt distributions and the characteristics of the hard and the soft samples, and it does not 



change the shape of the correlations. With the new Et threshold the fraction of K® per event rises 
by the same amount, around 30%, in the two samples. This means that the ratio of the rate of 
in soft events to the same rate in hard events is not influenced by the higher threshold. 



C. Analysis Discussion 

Some simple observations can be made about Table |1J The fraction of the total K® that falls 
into the soft subsample is rather small, ranging from about 30% at 630 GeV to about 18% at 1800 
GeV (19% and 10% for A respectively). The corrected mean number of K® produced per event 
in the full MB sample is about (8.6±2.6)% at 630 GeV and (8.8±2.6)% at 1800 GeV (respectively 
(3.7±1.1)% and (4.3±1.0)% for A ). 

The /it cross-section ratio may be obtained by fitting the K® and charged-track invariant pt 
distributions in the available pt range and extrapolating the fitted functions down to the minimum 
Pt value. This ratio is evaluated both for pTmin=0 GeV/ c and PTmin=0.4 GeV/ c. With the 
above technique, a ratio of 0.13±0.04 (including the systematic uncertainty) at -^5=1800 GeV and 
0.18±0.05 at 630 GeV is obtained for PTmin—OA GeV/ c. The same ratio for pTmin=0 GeV/ c, gives 
0.14±0.05 at v^=1800 GeV and 0.19±0.06 at 630 GeV (Table EJ- These last measurements are 
compatible with the previous CDF results Q], though are slightly higher at a/s = 630 GeV. In 
Table |HJ the corresponding values for the soft and hard subsamples are reported. It is remarkable 
that Kg I tt ratio is about two times larger in soft than in MB events. 

Studies of the production of strange particles K® and A in proton-antiproton interactions at 
different y/s are described in Ref. [ij at ^=540 GeV and |lj at 200 and 900 GeV. In Refs. |n|, 
(21 1 . results at a/s=1800 GeV are presented. Comparison with our results is restricted to the 



Q, 

full MB samples; furthermore, it should be noted that here no absolute cross- sect ions are provided. 
Comparison with Refs. ^^-[^ also requires taking into account the different pt and rj regions 
selected. 



A direct comparison of the invariant pt distribution of K® can be done with Ref. [17|. There the 
Pt distribution of K® is fitted to the functional power-law form of Eq.(0), fixing the parameter p° 
to 1.3 GeV/c. The average pr is computed from the parameters of the fit as: 

W = 2^- (4) 
n — 3 

With the new increased statistics and larger pt range, the fit with the p° parameter fixed, while 



giving a reasonable description of the px spectrum in the \ow-px region, does not describe the 
data at higher px- It yields a px compatible with the previous one (see Table EJ but with a large 
X 2 '/Nd.o.f.- The best fit to our distribution (shown in Figures 171 to fTOl for the soft and MB samples) 
is obtained with this form when all three parameters are allowed to vary freely and the fit region is 
restricted to px >1 GeV/ c. A summary of the results is reported in Table IV! The measurements 
reported in Refs. Ql-jsjj were done at different energies and in different phase space regions. 

From our best fit of MB sample data at 1800 GeV (630 GeV), the mean px of K® is 
(0.75 ± 0.07) GeV/ c ((0.70 ± 0.08) GeV/ c). These values are significantly higher than the previous 
CDF measurements due to the higher statistics in the high px tail of the distribution. 

Taking into account the different conditions and the method of measurement, it is possible to 

compare to UA5 data (Refs. [3], ^|) as well; our present measurement is also higher in this case. 

For completeness, the fit results of the px invariant distribution for the soft subsample are also 

reported in Table A second fitting function used is of the form: 

d?Nk 2 
E = exp(A + Bpx), where px = —— ■ (5) 
dp 6 B 

At both energies, we obtain a good x 2 /^d.o.f. using this function (see Table Ej). Therefore, the 
shape of the soft distribution is also well described by an exponential function; the mean px of the 
fit is generally larger than what is obtained using Eq.(J2J). For A , a systematically higher mean px 
than other experiments at equivalent energy is obtained (compare with Refs. and 20|). In this 
case as well, MB data can be equally well fitted by form (j2J) and by an exponential function. A 
summary of these results is in Table IVII 

The increase of the mean px (computed as in Eq. (j3J)) of the observed as a function of the 
event charged multiplicity is always larger than that of charged tracks. The increase for A is 
even larger, leading to the conclusion that it depends on the particle mass, as expected. A similar 



analysis is also reported in Ref. 



2l| . A direct comparison is not possible because of the different 



Px range and 77 acceptance, which reflect in larger multiplicities. However, a rise in mean px with 
heavier particle masses is clearly observed. 

In the analysis of charged tracks [4], all the correlations examined in the MB and in the hard 
samples showed different behaviors with respect to E cms , while a clear invariance was seen in the 
soft sample. With the available iff (A ) statistics it is not possible to discern any difference in 
the (px) dependence on multiplicity at the two energies, even in the full MB sample. Neverthe- 
less, the behavior of the three subsamples is clearly different. We note that the mean K®(A°) 



Pt increases with N*p also in the soft subsample, a feature that is not explained by the current 
Si P] Q Q !3) This observation also holds for charged hadrons, as discussed in j^j]. 



mo dels ([ 

The ratios of the mean numbers of K®(A°) per event to the charged multiplicity drop in the 
first few bins (0 < N* h < 6) and are roughly constant for N* h > 6 (MB sample) for both K® and 
A . The dependence on N* h is more pronounced for K® than for A . The fraction of A per event 
and per track is obviously smaller than that of K® and for both is larger at 630 GeV than at 1800 
GeV. Finally, the dependencies of the number of A for the soft and the MB samples on N* h , 
besides differing by about a factor of two, are both roughly flat and different in shape from the 
corresponding of the distributions. 



VIII. CONCLUSIONS 



The present measurements extend the studies of charged particle properties in MB pp interactions 
to and A production. Using the data available at the two c.m.s. energies obtained under 
the same experimental conditions and similar statistics, we are able to directly compare the V° 
production properties at the two c.m.s. energies. Our results offer new findings and significant 
improvements to the existing knowledge of V° production. We summarize our results as follows: 

• The overall production rates of K® and A are in agreement with previous measurements. 

• The inclusive pt spectra of K® and A now extend to pt > 8 GeV/ c. The K® distribution 
shows a more detailed shape in the high p? region when compared to previous data. For both 

and A , we measure an average px significantly higher than previous results. 

• New results are presented on the distribution of K® and A multiplicity. 

• For the first time, the MB sample has been used to analyze V° production properties in its 
soft and hard components. Inclusive pt and multiplicity distributions of V° are shown for 
the soft and hard data. 

• Analyses of the dependence of the mean ff s °(A°) pt with the event charged multiplicity are 
presented. The observed dependence is not explained by the current theoretical models. 
Comparison with an analogous study performed on charged tracks indicates that the rate of 
the dependence grows with particle mass. An increase of the mean pt is observable also in 
the soft subsample alone. 



• The event charged multiplicity has been adopted as the independent variable to analyze the 
ratio of the mean number of K®(A°) per event to the number of primary charged particles. 
For both K® and A this ratio rises toward very low multiplicity, remaining roughly constant 
for N* h > 5. 
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TABLE I: Raw and corrected numbers of K® and A found in each data set. In the rightmost three columns 
the fraction of Kg (A ) per event is shown. The uncertainties on all the corrected numbers and fraction of 
A" S °(A°) per event include the systematic uncertainty. The total number of MB events at ^=1800 GeV 
(630 GeV) is 2,079,558(1,963,157). 



RAW 

MB Soft Hard MB (xlO 3 ) 


CORRECTED 
Soft (xlO 3 ) Hard (xlO 3 ) 


FRACTION of ff°(A°)/EVENT (%) 
MB Soft Hard 


1800 K° 36642 6733 29909 180±50 
GeV A 7518 782 6736 90±30 


34±10 150±40 
9±3 80±20 


8.8±2.6 3.5±1.0 
4.3±1.0 1.0±0.3 


13.3±4.0 
7.2±2.2 


630 K® 32222 9835 22387 170±50 
GeV A 5883 1098 4785 70±20 


50±15 120±35 
13±4 60±20 


8.6±2.6 4.5±1.4 
3.7±1.1 1.2±0.4 


14.3±4.3 
7.1±2.1 


TABLE II: K® /it ratio in each data set. Data computed in the full px range and for px> 0.4 GeV/c are 
shown. The ratios are evaluated integrating the px distributions by extrapolating the fitted function down 
to PTmin=0- (0.4) GeV/ c. Here efficiency corrections and systematic uncertainties are included. 


PTmin (GeV/c) 


MB 


Soft 


Hard 


1800 0.0 

0.4 


0.14±0.05 
0.13 ± 0.04 


0.38±0.12 
0.30 ± 0.09 


0.11±0.04 
0.11 ± 0.03 


630 0.0 

0.4 


0.19±0.06 
0.18 ± 0.05 


0.42±0.13 
0.33 ± 0.10 


0.14±0.05 
0.18 ± 0.06 



TABLE III: Summary of all systematic uncertainties. For each systematic uncertainty source its effect 
on the various measured quantities is reported. The symbol "-" means no effect on the corresponding 
quantity. 



Systematic 
uncertainty 



source 



(p T ) vs Mult. PT distr. N(V°) distr. <N(F°)> vs Mult N(V°) K/ir 
(included (— — - not included in figures - — ) (included 

in figures) in table [H]) 



MC simulation 



3% 10% 
4% (A ) 



10-25% 



2-20% 



IjOW-PT 

extrapolation 



3-20 



5% 

10% (A ) 



5% (K°) 
10% (A ) 



5% (K°) 
10% (A ) 



Primary vertex 
selection 



< 9% 



Global normalization 
factor 



< 30% 



TABLE IV: Number of events with N K (N A A ) divided by the total number of events. MB, soft and 
hard data at 1800 and 630 GeV are reported. 



N K/A — 1800 GeV - — — 630 GeV 



MB (xl0~ 3 ) Soft (xlO^ 3 ) Hard(xl0" 3 ) MB (xl0~ 3 ) Soft (xl0~ 3 ) Hard (xl0~ 3 ) 








919± 7 


966± 2 


880± 10 


920± 10 


956± 5 


870± 20 




1 


80± 10 


33± 6 


110± 20 


80± 20 


40± 10 


120± 30 


K° 


2 


5± 1 


0.9± 0.2 


9± 2 


4± 1 


1.2± 0.3 


8± 2 




3 


0.6± 0.2 




1.0± 0.3 


0.3± 0.1 


0.04± 0.02 


0.7± 0.3 




4 


0.2± 0.1 




0.4± 0.2 













958± 5 


990.3± 0.1 


930± 10 


964± 5 


988.3± 0.1 


930± 10 


A 


1 


40± 20 


10± 4 


70± 30 


40± 20 


12± 5 


70± 30 




2 


1.6± 0.8 


0.08± 0.05 


3± 1 


0.9± 0.5 




2± 1 



TABLE V: Results of the fit to the invariant pt distribution of Kg. Data at different E cms are reported 
(for different experiments the parameters of the fit were reported when available). Parameters p° and n 
refer to the power-law (P.L.) function (Eq.(J2J)], B to the exponential (Exp.) form [Eq.l(5))]. CDF-0 refers 
to the so called Run-0 of the Tevatron and CDF-I to Run I data (this analysis). 



Experiment 


Data Set 


Pt 


p° (P.L.) 


n (P.L.) 


B (Exp.) 


X 2 /N d . oJ . 


(y/s in GeV) 




(GeV/c) 


(GeV/c) 








UA5(*) (546) [IS] 


MB 


0.58±0.04 


— 


— 




1.15 


CDF-0 (630) [17] 


MB 


0.5±0.1 


1.3(fixed) 


7.9±0.03 




3.9 


CDF-I (630) 


MB 


0.70±0.08 


3.3±0.2 


12.6±0.6 




68/57 


UA5(*) (900)[19j 


MB 


0.63±0.03 








0.5 


CDF-0 (1800) [17] 


MB 


0.60±0.03 


1.3(fixed) 


7.7±0.2 




0.74 


CDF-I (1800) 


MB 


0.58±0.02 


1.3(fixed) 


7.49±0.02 




265/68 


CDF-I (1800) 


MB 


0.75±0.07 


3.29±0.08 


11.7±0.1 




67/67 


CDF-I (630) 


Soft 


0.58±0.04 


9.0±0.1 


33.7±0.1 




29/22 


CDF-I (630) 


Soft 


0.64±0.02 






-3.12±0.03 


24/23 


CDF-I (1800) 


Soft 


0.62±0.02 


9.5±0.3 


33.7±0.9 




23/25 


CDF-I (1800) 


Soft 


0.67±0.02 






-3.00±0.04 


29/26 



(*) UA5 fits to a power-law form in px >0.4 together with an exponential form in px <0.4 GeV/c. 



TABLE VI: Results of the fit to the invariant pt distribution of A . Data at different E cms are reported 
(for different experiments the parameters of the fit were reported when available). Parameters p° and n 
refer to the power-law (P.L.) function [Eq.@], B to the exponential (Exp.) form [Eq.l|5|)]. CDF-I refers to 
Run I data (this analysis). 



Experiment 
(y/s in GeV) 


Data Set 


Pt 
(GeV/c) 


p° (P.L.) 

(GeV/c) 


n (P.L.) B (Exp.) 


X 2 /N d . . f . 


UA5 (546) 


MB 


0.62±0.08 








CDF-I (630) 


MB 


0.91±0.07 


12.3±0.1 


30.1±0.2 


59/35 


CDF-I (630) 


MB 


0.98±0.01 




-2.05±0.03 


50/36 


UA5 (900) [12] 


MB 


0.97±0.01 








CDF-I (1800) 


MB 


0.97±0.09 


12.4±0.1 


28.6 ± 0.09 


41/45 


CDF-I (1800) 


MB 


1.04±0.01 




-1.92±0.02 


55/46 


CDF-I (630) 


Soft 


0.67±0.09 


10.0±0.2 


33.0±0.2 


31/24 


CDF-I (630) 


Soft 


0.73±0.1 




-2.74±0.05 


28/25 


CDF-I (1800) 


Soft 


0.64±0.05 


9.5±3.3 


33.0 ±0.2 


29/22 


CDF-I (1800) 


Soft 


0.73±0.10 




-2.74±0.05 


25/23 
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FIG. 1: Invariant mass with a ir — ir and 7r-p mass assignment for oppositely signed track pairs passing all 
selection requirements but allowing the mass to be within 10a from the nominal mass. The background 
includes the contamination by K® in the A sample and vice versa. 1800 GeV data are shown. 




FIG. 2: Overall correction factors at 1800 GeV for K°(A°) function of the transverse momentum. The 
correction factors are defined in the text. 




FIG. 3: Raw and corrected K® and A proper lifetime distributions for 1800 GeV data. The line represents 
the best exponential fit to the data. 
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FIG. 4: The mean j»t of K® and A as a function of proper lifetime (t) at 1800 GeV. Raw and corrected 
data are shown. 
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FIG. 5: Distribution of the multiplicity of Kg at 1800 (full symbols) and 630 GeV (open symbols). P(Nk) 
= (number of events with Nk i^ s °)/(total number of events). MB, soft (divided by 100) and hard data 
(multiplied by 100) are shown. 
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FIG. 6: Distribution of the multiplicity of A at 1800 (full symbols) and 630 GeV (open symbols). P(N\) 
= (number of events with N\ A ) /(total number of events). MB, soft (divided by 100) and hard data 
(multiplied by 100) are shown. 
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FIG. 7: Kg inclusive invariant pt distributions at 1800 GeV. MB, soft and hard data are shown, normalized 
to the number of events in each sample. E is the particle energy and N even t is the total number of events 
which contribute to the distribution. To separate the curves, hard data points are multiplied by 10 and 
soft data points by 0.1. The solid lines represent the best fits to Eq.© of the text. 
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FIG. 8: Kg inclusive invariant pt distributions at 630 GeV MB, soft and hard data are shown, normalized 
to the number of events in each sample. E is the particle energy and N even t is the total number of events 
which contribute to the distribution. To separate the curves, hard data points are multiplied by 10 and 
soft data points by 0.1. The solid lines represent the best fits to Eq.© of the text. 
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FIG. 9: A inclusive invariant px distributions at 1800 GeV. MB, soft and hard data are shown, normalized 
to the number of events in each sample. E is the particle energy and N even t is the total number of events 
which contribute to the distribution. To separate the curves, hard data points are multiplied by 10 and 
soft data points by 0.1. The solid lines represent the best fits to Eq.© of the text. 
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FIG. 10: A inclusive invariant px distributions at 630 GeV. MB, soft and hard data are shown, normalized 
to the number of events in each sample. E is the particle energy and N even t is the total number of events 
which contribute to the distribution. To separate the curves, hard data points are multiplied by 10 and 
soft data points by 0.1. The solid lines represent the best fits to Eq.© of the text. 
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FIG. 11: Average transverse momentum (px) of K® and A at 1800 GeV as a function of the event charged 
multiplicity (N* h ). MB data are shown. For comparison, the mean pt of charged particles measured in 
the same phase space region is also plotted P| . The filled squares around the points delimit the systematic 
uncertainties. 
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FIG. 12: Average transverse momentum (px) of i^ s ° and A at 1800 GeV as a function of the event charged 
multiplicity {N* h ). Soft data are shown. For comparison, the mean pt of charged particles measured in 
the same phase space region is also plotted ^ | . The filled squares around the points delimit the systematic 
uncertainties. 



1.8 \-n A 
• K' 



.6 



> 

<J> 1 .4 
O 



A 

H 
Oh 

V 



1 .2 

1 

0.8 
0.6 



Hard 
o 



800 GeV 



o 



tracks 



I* 



V 



[771^1 .0 

p T ^0.4 GeV/c 



5 10 15 20 25 30 35 40 



FIG. 13: Average transverse momentum (px) of i^ s ° and A at 1800 GeV as a function of the event charged 



multiplicity (N* h ). Hard data are shown, 
the same phase space region is also plotted 
uncertainties. 



or comparison, the mean pj- of charged particles measured in 
4] . The filled squares around the points delimit the systematic 




FIG. 14: Average transverse momentum (jot) of and A at 630 GeV as a function of the event charged 
multiplicity (N* h ). MB data are shown. For comparison, the mean pt of charged particles measured in 
the same phase space region is also plotted P| . The filled squares around the points delimit the systematic 
uncertainties. 
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FIG. 15: Average transverse momentum (jot) of and A at 630 GeV as a function of the event charged 
multiplicity {N* h ). Soft data are shown. For comparison, the mean pt of charged particles measured in 
the same phase space region is also plotted ^ | . The filled squares around the points delimit the systematic 
uncertainties. 
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FIG. 16: Average transverse momentum (jot) of and A at 630 GeV as a function of the event charged 



multiplicity (N* h ). Hard data are shown, 
the same phase space region is also plotted 
uncertainties. 



^or comparison, the mean pj- of charged particles measured in 
4] . The filled squares around the points delimit the systematic 
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FIG. 17: Mean number of K? per event divided by the charged multiplicity (N* h ) and plotted as a function 
of N* h . The MB, so/t and hard data at 1800 GeV are shown. 
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FIG. 18: Mean number of K° per event divided by the charged multiplicity {^*h) an d plotted as a 
function of N* h . The MB, soft and /iarcf data at 1800 GeV are shown. 
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FIG. 19: Mean number of A per event divided by the charged multiplicity (N* h ) and plotted as a function 
of N* h . The MB, soft and hard data at 1800 GeV are shown. 
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FIG. 20: Mean number of A per event divided by the charged multiplicity (N* h ) and plotted as a function 
of N* h . The MB, soft and hard data at 630 GeV are shown. 



